Entropic transport: Kinetics, scaling and control mechanisms 
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We show that transport in the presence of entropic barriers exhibits peculiar characteristics which 
makes it distinctly different from that occurring through energy barriers. The constrained dynamics 
yields a scaling regime for the particle current and the diffusion coefficient in terms of the ratio 
between the work done to the particles and available thermal energy. This interesting property, 
genuine to the entropic nature of the barriers, can be utilized to effectively control transport through 
quasi one-dimensional structures in which irregularities or tortuosity of the boundaries cause entropic 
effects. The accuracy of the kinetic description has been corroborated by simulations. Applications 
to different dynamic situations involving entropic barriers are outlined. 
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Transport through quasi-onedimensional structures as 
pores, ion channels and zeolites is ubiquitous in biologi- 
cal and physico-chemical systems and constitute a basic 
mechanism in processes as catalysis, osmosis and particle 
separation 0, 0, H El El • A common characteristic of 
these systems is the confinement arising from the pres- 
ence of boundaries which very often exhibit an irregular 
geometry. Variations of the shape of the structure along 
the propagation direction implies changes in the number 
of accessible states of the particles. Consequently, en- 
tropy is spatially varying, and the system evolves through 
entropic barriers, which controls the transport, promot- 
ing or hampering the transfer of mass and energy to cer- 
tain regions. Motion in the system can be induced by 
the presence of external driving forces supplying the par- 
ticles with the energy necessary to proceed. The study of 
the kinetics of the entropic transport, the properties of 
transport coefficients in far from equilibrium situations 
and the possibility for transport control mechanisms are 
objectives of major importance in the dynamical charac- 
terization of those systems. 

Our purpose in this Letter is to demonstrate that 
entropic transport exhibits striking features, sometimes 
counterintuitive, which are different from those ob- 
served in the more familiar case with energy barriers • 
We propose a general scenario describing the dynamics 
through entropic barriers and show the existence of a 
scaling regime for the current of particles and the effective 
diffusion coefficient. The presence of this regime might 
have important implications in the control of transport. 

Entropic transport. - The origin of the entropic barri- 
ers can be inherent to the intimate nature of the system 
or may emerge as a consequence of a coarsening of the 
description employed. A typical example presents the 
motion of a Brownian particle in an enclosure of varying 
cross-section. This basic situation constitutes the start- 
ing point in the study of transport processes in the type 
of confined systems that are very often encountered at 
sub-cellular level, nanoporous materials and in microflu- 
idic applications. As shown in Ref. the compli- 
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FIG. 1: Schematic diagram of the tube confining the motion of 
the biased Brownian particles. The half width lo is a periodic 
function of x with periodicity L. 

cated boundary conditions of the diffusion equation in 
irregular channels can be greatly simplified by introduc- 
ing an entropic potential that accounts for the reduced 
space accessible for the diffusion of the Brownian particle. 
The resulting kinetic equation describing the evolution of 
the probability distribution, is known as the Fick- Jacobs 
equation [8|, |9j and constitutes an approximation to the 
full dynamics. The validity of that equation has only 
been analyzed for diffusion in the absence of a driving 
force in whose case many of the transport processes pre- 
viously mentioned could not take place. This is so, since 
thermal diffusion alone may not be able to induce tran- 
sitions of the particles through the entropic barrier. 

In typical transport processes through pores or chan- 
nels, motion of the suspended particles is induced by ap- 
plication of an external driving force F that is directed 
along their axis. The over-damped dynamics of a biased 
Brownian particle within the tube (see Fig.^l then reads: 

rj— =F + v^feTf(*), (1) 

where rj is the friction coefficient of the particle, fee the 
Boltzmann constant, T the temperature, F — \F\ a con- 
stant force in x-direction and £(t) is Gaussian white noise 
with zero mean and correlation function: (£i(t) = 
2 6ij6(t — t') for i,j — x,y,z. The reflecting boundary 
conditions ensure the confinement of the dynamics within 
the tube. 
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Reduction of the dimensionality. - As mentioned pre- 
viously, the dynamics of the particles along the axis of 
the 3D tube or a 2D channel (see Fig. Q can be recast 
into the Fick- Jacobs equation; i.e., 
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obtained from the 3D (or 2D) Smoluchowski equation af- 
ter elimination of y and z coordinates by assuming equi- 
librium in the orthogonal directions. Here P(x,t) is the 
probability distribution function, Dq the diffusion coef- 
ficient, and s(x) is the cross-sectional area for a (3D) 
tube or the width for a (2D) channel. This description 
is in principle valid for |o/(a;)| <C 1, where u>(x) is the 
radius of the tube (or the halfwidth of the channel in 
2D) and the prime refers to the first derivative. It has 
been shown that the introduction of a ^-dependent diffu- 
sion coefficient considerably improves the accuracy of the 
kinetic equation extending its validity to more winding 
structures @, El ■ The expression 
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where Do — ksT/r] and a — 1/3, 1/2 for two and three 
dimensions, respectively, was shown to appropriately ac- 
count for the curvature effects [Icj . 

In the presence of a constant force F along the direc- 
tion of the tube the Fick- Jacobs equation can be recast 
into the following expression ^3 
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which defines a free energy A(x) := E — T S = —Fx — 
Tfce hih(x), where E = —Fx is the energy, S — 
k-Q hih(x) the entropy, h{x) the dimensionless width 
2w(x)/L in 2D, and the dimensionless transverse cross- 
section 7r [uj(x)/L] 2 of the tube in 3D. For a symmetric 
channel with periodicity L, the free energy assumes the 
form of a periodic tilted potential. 

Universal scaling for the particle current and effective 
diffusion. - The key quantities in transport through 
quasi one-dimensional structures are the average particle 
current, or equivalently the nonlinear mobility, and the 
effective diffusion coefficient. While in the case of an 
energy barrier, the driving force F and the temperature 
T are two independent variables, for entropic transport, 
both current and effective diffusion are controlled by a 
universal scaling parameter: 
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For the average particle current and the nonlinear mobil- 
ity n(f) we find an expression similar to the Stratonovich 
formula 
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where 



/(*,/):= l + ^'(z)f h~\z) 



x exp(-/z) / dz h(z) exp(f z) , (6b) 

Jz-l 

depends only on the dimensionless variable z = x/L, the 
scaling parameter / and the shape of the tube given in 
terms of the dimensionless half width uj(z) := w(x)/L 
and its first derivative. Here h(z) := h{x). 

The effective diffusion coefficient could be expressed in 
terms of moments of the first passage time for a Brownian 
particle arriving at xq + L while starting out from xo 
[ll| . A detailed analysis shows that the effective diffusion 
coefficient also scales with F ' Ljk^T as: 
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Numerical simulations. - A model of a 2D periodic 
channel is sketched in Fig. ^ the shape is described by 
lu(x) = a sm(2irx/ L) + b. Here, a is the parameter that 
controls the slope of the walls, the width of the channel 
is 2lu(x), and the width at the bottleneck is 2 (b — a). 

The scaling behaviors, predicted above, have been cor- 
roborated by Brownian dynamic simulations performed 
by integration of the dimensionless Langevin equation, 
which is equivalent to Eq. lfT|l. within the stochastic 
Euler-algorithm. Therefore lengths are scaled by the pe- 
riodicity L of the tube, time by r := L 2 ry/(/cB T room ) - 
the corresponding characteristic diffusion time at room 
temperature T Toom - and the force by F := ijL/t. The 
mean velocity in ^-direction, (x) — limt—,.^ x(t)/t, and 
the corresponding effective diffusion coefficient, D e g = 
1/2 lim t ^ 00 ((a; 2 (t)) — (x(t)) 2 )/t, are obtained as an av- 
erage over 3 • 10 4 trajectories. 

Results for the particle current and the effective dif- 
fusion coefficient as a function of the applied force for 
the case a — 1/(2 7r), b = 1.02/(2 if) and L — 1 are pre- 
sented in Fig.|2|and Fig.[3]for different values of the noise 
strength (i.e. the temperature). The particle current in- 
creases monotonically with the force, but decreases upon 
increasing the level of noise. The effective diffusion coeffi- 
cient exhibits a non-monotonic behavior with the appear- 
ance of a peak which becomes more pronounced at low 
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FIG. 2: Numerically determined force dependence of particle 
current for a symmetric two-dimensional tube with the shape 
defined by the half width lj(x) = [sin(27ra;/L) + 1.02]/(2tt), 
L = 1 and for the values of T/T room : 0.01 (solid line), 0.1 
(dashed line), 0.2 (dotted line) and 0.4 (dash-dotted line). 
The inset depicts the dependence of the particle current 
(x)/(L/r) on the dimensionless temperature T/Troom for the 
force value: F/Fo = 0.628. Contrary to the case of energetic 
barriers, the particle current declines with increasing temper- 
ature. 



FIG. 4: Graph for the scaled nonlinear mobility. In the 
Langevin simulation the different symbols correspond to dif- 
ferent values of T/T room : 0.01 (crosses), 0.1 (pluses), 0.2 
(squares), 0.4 (triangles). The relative error of the simulation 
results is 0.01. The Fick-Jacobs results, Eq. ©, correspond 
to the solid lines. The inset depicts the long range behavior 
(the dotted line depicts the numerical results) . The numerical 
values for the scaled nonlinear mobility approach, in the limit 
/ — > oo, to the value 1 (dashed horizontal line). 




FIG. 3: The effective diffusion coefficient versus the exter- 
nal bias. The parameters for the various lines correspond to 
those detailed in Fig. |2] The inset depicts the effective dif- 
fusion coefficient D e g/(L 2 /r) vs. dimensionless temperature 
T/T TO 



noise levels (see Fig. EJ . When both quantities are repre- 
sented as a function of the scaling parameter / (see Fig.0] 
and Fig.[5j) all curves collapse to the scaled solution which 
evidences the excellent agreement of simulations results 
with the scaling behavior predicted for those quantities. 
Therefore, whereas in the case of transport through en- 
ergetic barriers the force (or tilt) and the temperature 
are two independent parameters, the entropic transport 
is controlled by a single parameter /. Another important 
results shown in Fig. [3] is the presence of a peak in the 
diffusion and the fact that the effective diffusion can be 
much larger than bulk diffusion. Thus the phenomenon 
of enhancement of the diffusion, linked to the dynamics 
of particles in periodic tilted energetic potentials, also 



takes place when barriers hindering the transport have 
an entropic nature. 

In Fig. 0] and Fig. [5] we have also represented the 
nonlinear mobility and effective diffusion coefficient pre- 
dicted by Eq. ijfJJl and (JJJ obtained from the Fick-Jacobs 
equation. At low values of the scaling variable / the re- 
sults match perfectly with the simulations whereas devi- 
ations occur at higher values of /. The scaled nonlinear 
mobility and the effective diffusion coefficient approxi- 
mates for / — > oo values different from the value 1. The 
accuracy of the Fick-Jacobs description worsens at large 
/ because the assumption of equilibration in the trans- 
verse direction, which supports the elimination of the y,z 
coordinates, fails at high values of the applied force. The 
agreement substantially improves when the shape of the 
tube does not change too fast, i.e. when |u/(ir)| is smaller, 
which can be achieved for instance by increasing the pe- 
riod L of the shape oscillations of the channel. In situa- 
tions where the roughness of the channel is not very ex- 
treme, the Fick-Jacobs description provides a very good 
approximation to the transport for values of the external 
work of some tens of fee T's. In fact, that is the range of 
energies relevant to most transport processes in biological 
systems. 

The peculiar behavior of nonlinear mobility and effec- 
tive diffusion coefficient as a function of temperature is 
depicted in the insets of Fig. [5] and [3] Contrary to the 
case of an energetic barrier, the particle current decreases 
upon increasing the temperature. In the presence of ener- 
getic barriers, the temperature facilitates the activation 
(the overcoming of the barriers) and thus tends to in- 
crease the particle current. However, when transport is 
controlled by entropic factors, the temperature dictates 
the strength of the entropic potential, and thus an in- 
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FIG. 5: Same as in Fig.0] but for the effective diffusion coef- 
ficient. The relative error of the simulation results is 0.1. The 
Fick- Jacobs results, Eq. (J, correspond to the solid lines. 
The inset depicts the long range behavior (the dotted line 
depicts the numerical results). The numerical values for the 
scaled nonlinear mobility approach, in the limit / — > oo, to 
the value 1 (dashed horizontal line) . 



crease of temperature leads to a reduction of the particle 
current. The effective diffusion coefficient as a function 
of the temperature also manifest a striking behavior with 
the presence of a peak, and the existence of a range of 
temperatures where the effective diffusion coefficient de- 
creases upon increasing the temperature. It is important 
to remark that, since the transport characteristics scale 
as FL/kBT, the peculiar regimes can be obtained not 
only by changing the temperature but also by modifying 
the strength of the force. 

Applications. - An example in which the entropic na- 
ture of the transport becomes more evident is the case of 
micro and nanoporous materials, such as zeolites. These 
materials have a regular structure with channels of dif- 
ferent width and well-defined geometry. This peculiar 
structure confers them an outstanding ability to act as 
molecular sieves, that is currently exploited in chemically 
clean separation of mixtures, ion exchange and petro- 
chemical cracking. Driven by their economic and scien- 
tific importance, these materials have been studied ex- 
tensively experimentally and more recently by computer 
simulations. For instance, the diffusion has been found to 



decrease with temperature in some range of temperatures 
fl2| ; and the existence of an optimal value of the diffusion 
as a function of the temperature has also been observed 
fl3|. In fact, the dependence of the effective diffusion 
coefficient on temperature reported in Ref. behaves 
just as the one predicted here with Fig.OH Finally, values 
of diffusion coefficients higher than the bulk, consistent 
with the phenomenon of diffusion enhancement predicted 
by our model, have also been reported Q. Our simple 
model thus accounts for all these behaviors and shows 
that they are not specific of a particular zeolite structure 
but they arise from the entropic nature of the transport. 

Conclusions. - In summary, we have shown that 
transport phenomena in systems in the presence of en- 
tropic barriers exhibit some features radically different 
from conventional transport through energetic barriers. 
The effect of confinement can be recast in terms of an en- 
tropic potential, and the dynamics of the system can be 
accurately described by means of the Fick- Jacobs equa- 
tion. Wc have shown the existence of a scaling regime 
in the dynamics. The particle current and the effective 
diffusion coefficient are controlled by a single parameter 
/ that measures the relative importance of the exter- 
nal work done to the particle and the thermal energy. 
The scaling in / thus opens up the possibility of tuning 
and controlling the efficiency of transport in confined sys- 
tems by a proper combination of temperature and applied 
field. In situations in which the temperature can only be 
varied in a very limited range, as frequently occurs in 
biological systems, the existence of scaling implies that 
the same transport regime can be accomplished by the 
application of an external force. The analysis presented 
could be applied to a wide variety of situations, such as 
biological transport through ion channels and membrane 
pores, or the portage in molecular sieves or polymer gels, 
where entropic effects play a very important role. 
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